\(\int \frac {1}{x^4 (a^2+2 a b x^3+b^2 x^6)^{3/2}} \, dx\) [106]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (verified)
   Maple [C] (warning: unable to verify)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [A] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [F(-1)]

Optimal result

Integrand size = 26, antiderivative size = 188 \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=-\frac {2 b}{3 a^3 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {b}{6 a^2 \left (a+b x^3\right ) \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {a+b x^3}{3 a^3 x^3 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {3 b \left (a+b x^3\right ) \log (x)}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}}+\frac {b \left (a+b x^3\right ) \log \left (a+b x^3\right )}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}} \]

[Out]

-2/3*b/a^3/((b*x^3+a)^2)^(1/2)-1/6*b/a^2/(b*x^3+a)/((b*x^3+a)^2)^(1/2)+1/3*(-b*x^3-a)/a^3/x^3/((b*x^3+a)^2)^(1
/2)-3*b*(b*x^3+a)*ln(x)/a^4/((b*x^3+a)^2)^(1/2)+b*(b*x^3+a)*ln(b*x^3+a)/a^4/((b*x^3+a)^2)^(1/2)

Rubi [A] (verified)

Time = 0.07 (sec) , antiderivative size = 188, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.115, Rules used = {1369, 272, 46} \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=-\frac {b}{6 a^2 \left (a+b x^3\right ) \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {3 b \log (x) \left (a+b x^3\right )}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}}+\frac {b \left (a+b x^3\right ) \log \left (a+b x^3\right )}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {2 b}{3 a^3 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {a+b x^3}{3 a^3 x^3 \sqrt {a^2+2 a b x^3+b^2 x^6}} \]

[In]

Int[1/(x^4*(a^2 + 2*a*b*x^3 + b^2*x^6)^(3/2)),x]

[Out]

(-2*b)/(3*a^3*Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6]) - b/(6*a^2*(a + b*x^3)*Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6]) - (a +
b*x^3)/(3*a^3*x^3*Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6]) - (3*b*(a + b*x^3)*Log[x])/(a^4*Sqrt[a^2 + 2*a*b*x^3 + b^2*
x^6]) + (b*(a + b*x^3)*Log[a + b*x^3])/(a^4*Sqrt[a^2 + 2*a*b*x^3 + b^2*x^6])

Rule 46

Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d*x
)^n, x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - a*d, 0] && ILtQ[m, 0] && IntegerQ[n] &&  !(IGtQ[n, 0] && Lt
Q[m + n + 2, 0])

Rule 272

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a
+ b*x)^p, x], x, x^n], x] /; FreeQ[{a, b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]

Rule 1369

Int[((d_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_.) + (c_.)*(x_)^(n2_.))^(p_), x_Symbol] :> Dist[(a + b*x^n + c*x^
(2*n))^FracPart[p]/(c^IntPart[p]*(b/2 + c*x^n)^(2*FracPart[p])), Int[(d*x)^m*(b/2 + c*x^n)^(2*p), x], x] /; Fr
eeQ[{a, b, c, d, m, n, p}, x] && EqQ[n2, 2*n] && EqQ[b^2 - 4*a*c, 0] && IntegerQ[p - 1/2]

Rubi steps \begin{align*} \text {integral}& = \frac {\left (b^2 \left (a b+b^2 x^3\right )\right ) \int \frac {1}{x^4 \left (a b+b^2 x^3\right )^3} \, dx}{\sqrt {a^2+2 a b x^3+b^2 x^6}} \\ & = \frac {\left (b^2 \left (a b+b^2 x^3\right )\right ) \text {Subst}\left (\int \frac {1}{x^2 \left (a b+b^2 x\right )^3} \, dx,x,x^3\right )}{3 \sqrt {a^2+2 a b x^3+b^2 x^6}} \\ & = \frac {\left (b^2 \left (a b+b^2 x^3\right )\right ) \text {Subst}\left (\int \left (\frac {1}{a^3 b^3 x^2}-\frac {3}{a^4 b^2 x}+\frac {1}{a^2 b (a+b x)^3}+\frac {2}{a^3 b (a+b x)^2}+\frac {3}{a^4 b (a+b x)}\right ) \, dx,x,x^3\right )}{3 \sqrt {a^2+2 a b x^3+b^2 x^6}} \\ & = -\frac {2 b}{3 a^3 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {b}{6 a^2 \left (a+b x^3\right ) \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {a+b x^3}{3 a^3 x^3 \sqrt {a^2+2 a b x^3+b^2 x^6}}-\frac {3 b \left (a+b x^3\right ) \log (x)}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}}+\frac {b \left (a+b x^3\right ) \log \left (a+b x^3\right )}{a^4 \sqrt {a^2+2 a b x^3+b^2 x^6}} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(901\) vs. \(2(188)=376\).

Time = 1.33 (sec) , antiderivative size = 901, normalized size of antiderivative = 4.79 \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=-\frac {-2 a^6-5 a^5 b x^3+2 a^4 b^2 x^6+4 a^3 b^3 x^9-a b^5 x^{15}+2 a^4 \sqrt {a^2} \sqrt {\left (a+b x^3\right )^2}+3 a^3 \sqrt {a^2} b x^3 \sqrt {\left (a+b x^3\right )^2}-5 \left (a^2\right )^{3/2} b^2 x^6 \sqrt {\left (a+b x^3\right )^2}+a \sqrt {a^2} b^3 x^9 \sqrt {\left (a+b x^3\right )^2}-\sqrt {a^2} b^4 x^{12} \sqrt {\left (a+b x^3\right )^2}+6 b x^3 \left (\left (a^2\right )^{3/2} b^2 x^6+a^4 \left (\sqrt {a^2}-\sqrt {\left (a+b x^3\right )^2}\right )+a^3 b x^3 \left (2 \sqrt {a^2}-\sqrt {\left (a+b x^3\right )^2}\right )\right ) \text {arctanh}\left (\frac {b x^3}{\sqrt {a^2}-\sqrt {\left (a+b x^3\right )^2}}\right )-6 b x^3 \left (a^5+2 a^4 b x^3-\left (a^2\right )^{3/2} b x^3 \sqrt {\left (a+b x^3\right )^2}+a^3 \left (b^2 x^6-\sqrt {a^2} \sqrt {\left (a+b x^3\right )^2}\right )\right ) \log \left (x^3\right )+3 a^5 b x^3 \log \left (\sqrt {a^2}-b x^3-\sqrt {\left (a+b x^3\right )^2}\right )+6 a^4 b^2 x^6 \log \left (\sqrt {a^2}-b x^3-\sqrt {\left (a+b x^3\right )^2}\right )+3 a^3 b^3 x^9 \log \left (\sqrt {a^2}-b x^3-\sqrt {\left (a+b x^3\right )^2}\right )-3 a^3 \sqrt {a^2} b x^3 \sqrt {\left (a+b x^3\right )^2} \log \left (\sqrt {a^2}-b x^3-\sqrt {\left (a+b x^3\right )^2}\right )-3 \left (a^2\right )^{3/2} b^2 x^6 \sqrt {\left (a+b x^3\right )^2} \log \left (\sqrt {a^2}-b x^3-\sqrt {\left (a+b x^3\right )^2}\right )+3 a^5 b x^3 \log \left (\sqrt {a^2}+b x^3-\sqrt {\left (a+b x^3\right )^2}\right )+6 a^4 b^2 x^6 \log \left (\sqrt {a^2}+b x^3-\sqrt {\left (a+b x^3\right )^2}\right )+3 a^3 b^3 x^9 \log \left (\sqrt {a^2}+b x^3-\sqrt {\left (a+b x^3\right )^2}\right )-3 a^3 \sqrt {a^2} b x^3 \sqrt {\left (a+b x^3\right )^2} \log \left (\sqrt {a^2}+b x^3-\sqrt {\left (a+b x^3\right )^2}\right )-3 \left (a^2\right )^{3/2} b^2 x^6 \sqrt {\left (a+b x^3\right )^2} \log \left (\sqrt {a^2}+b x^3-\sqrt {\left (a+b x^3\right )^2}\right )}{3 a^4 \sqrt {a^2} x^3 \left (a^2+a b x^3-\sqrt {a^2} \sqrt {\left (a+b x^3\right )^2}\right )^2} \]

[In]

Integrate[1/(x^4*(a^2 + 2*a*b*x^3 + b^2*x^6)^(3/2)),x]

[Out]

-1/3*(-2*a^6 - 5*a^5*b*x^3 + 2*a^4*b^2*x^6 + 4*a^3*b^3*x^9 - a*b^5*x^15 + 2*a^4*Sqrt[a^2]*Sqrt[(a + b*x^3)^2]
+ 3*a^3*Sqrt[a^2]*b*x^3*Sqrt[(a + b*x^3)^2] - 5*(a^2)^(3/2)*b^2*x^6*Sqrt[(a + b*x^3)^2] + a*Sqrt[a^2]*b^3*x^9*
Sqrt[(a + b*x^3)^2] - Sqrt[a^2]*b^4*x^12*Sqrt[(a + b*x^3)^2] + 6*b*x^3*((a^2)^(3/2)*b^2*x^6 + a^4*(Sqrt[a^2] -
 Sqrt[(a + b*x^3)^2]) + a^3*b*x^3*(2*Sqrt[a^2] - Sqrt[(a + b*x^3)^2]))*ArcTanh[(b*x^3)/(Sqrt[a^2] - Sqrt[(a +
b*x^3)^2])] - 6*b*x^3*(a^5 + 2*a^4*b*x^3 - (a^2)^(3/2)*b*x^3*Sqrt[(a + b*x^3)^2] + a^3*(b^2*x^6 - Sqrt[a^2]*Sq
rt[(a + b*x^3)^2]))*Log[x^3] + 3*a^5*b*x^3*Log[Sqrt[a^2] - b*x^3 - Sqrt[(a + b*x^3)^2]] + 6*a^4*b^2*x^6*Log[Sq
rt[a^2] - b*x^3 - Sqrt[(a + b*x^3)^2]] + 3*a^3*b^3*x^9*Log[Sqrt[a^2] - b*x^3 - Sqrt[(a + b*x^3)^2]] - 3*a^3*Sq
rt[a^2]*b*x^3*Sqrt[(a + b*x^3)^2]*Log[Sqrt[a^2] - b*x^3 - Sqrt[(a + b*x^3)^2]] - 3*(a^2)^(3/2)*b^2*x^6*Sqrt[(a
 + b*x^3)^2]*Log[Sqrt[a^2] - b*x^3 - Sqrt[(a + b*x^3)^2]] + 3*a^5*b*x^3*Log[Sqrt[a^2] + b*x^3 - Sqrt[(a + b*x^
3)^2]] + 6*a^4*b^2*x^6*Log[Sqrt[a^2] + b*x^3 - Sqrt[(a + b*x^3)^2]] + 3*a^3*b^3*x^9*Log[Sqrt[a^2] + b*x^3 - Sq
rt[(a + b*x^3)^2]] - 3*a^3*Sqrt[a^2]*b*x^3*Sqrt[(a + b*x^3)^2]*Log[Sqrt[a^2] + b*x^3 - Sqrt[(a + b*x^3)^2]] -
3*(a^2)^(3/2)*b^2*x^6*Sqrt[(a + b*x^3)^2]*Log[Sqrt[a^2] + b*x^3 - Sqrt[(a + b*x^3)^2]])/(a^4*Sqrt[a^2]*x^3*(a^
2 + a*b*x^3 - Sqrt[a^2]*Sqrt[(a + b*x^3)^2])^2)

Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.12 (sec) , antiderivative size = 92, normalized size of antiderivative = 0.49

method result size
pseudoelliptic \(-\frac {\left (-3 b \,x^{3} \left (b \,x^{3}+a \right )^{2} \ln \left (b \,x^{3}+a \right )+3 b \,x^{3} \left (b \,x^{3}+a \right )^{2} \ln \left (b \,x^{3}\right )+a \left (3 b^{2} x^{6}+\frac {9}{2} a b \,x^{3}+a^{2}\right )\right ) \operatorname {csgn}\left (b \,x^{3}+a \right )}{3 \left (b \,x^{3}+a \right )^{2} a^{4} x^{3}}\) \(92\)
risch \(\frac {\sqrt {\left (b \,x^{3}+a \right )^{2}}\, \left (-\frac {b^{2} x^{6}}{a^{3}}-\frac {3 b \,x^{3}}{2 a^{2}}-\frac {1}{3 a}\right )}{\left (b \,x^{3}+a \right )^{3} x^{3}}-\frac {3 \sqrt {\left (b \,x^{3}+a \right )^{2}}\, b \ln \left (x \right )}{\left (b \,x^{3}+a \right ) a^{4}}+\frac {\sqrt {\left (b \,x^{3}+a \right )^{2}}\, b \ln \left (-b \,x^{3}-a \right )}{\left (b \,x^{3}+a \right ) a^{4}}\) \(116\)
default \(-\frac {\left (18 b^{3} \ln \left (x \right ) x^{9}-6 \ln \left (b \,x^{3}+a \right ) b^{3} x^{9}+36 b^{2} a \ln \left (x \right ) x^{6}-12 \ln \left (b \,x^{3}+a \right ) a \,b^{2} x^{6}+6 b^{2} x^{6} a +18 a^{2} b \ln \left (x \right ) x^{3}-6 \ln \left (b \,x^{3}+a \right ) a^{2} b \,x^{3}+9 a^{2} b \,x^{3}+2 a^{3}\right ) \left (b \,x^{3}+a \right )}{6 x^{3} a^{4} {\left (\left (b \,x^{3}+a \right )^{2}\right )}^{\frac {3}{2}}}\) \(133\)

[In]

int(1/x^4/(b^2*x^6+2*a*b*x^3+a^2)^(3/2),x,method=_RETURNVERBOSE)

[Out]

-1/3*(-3*b*x^3*(b*x^3+a)^2*ln(b*x^3+a)+3*b*x^3*(b*x^3+a)^2*ln(b*x^3)+a*(3*b^2*x^6+9/2*a*b*x^3+a^2))*csgn(b*x^3
+a)/(b*x^3+a)^2/a^4/x^3

Fricas [A] (verification not implemented)

none

Time = 0.25 (sec) , antiderivative size = 119, normalized size of antiderivative = 0.63 \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=-\frac {6 \, a b^{2} x^{6} + 9 \, a^{2} b x^{3} + 2 \, a^{3} - 6 \, {\left (b^{3} x^{9} + 2 \, a b^{2} x^{6} + a^{2} b x^{3}\right )} \log \left (b x^{3} + a\right ) + 18 \, {\left (b^{3} x^{9} + 2 \, a b^{2} x^{6} + a^{2} b x^{3}\right )} \log \left (x\right )}{6 \, {\left (a^{4} b^{2} x^{9} + 2 \, a^{5} b x^{6} + a^{6} x^{3}\right )}} \]

[In]

integrate(1/x^4/(b^2*x^6+2*a*b*x^3+a^2)^(3/2),x, algorithm="fricas")

[Out]

-1/6*(6*a*b^2*x^6 + 9*a^2*b*x^3 + 2*a^3 - 6*(b^3*x^9 + 2*a*b^2*x^6 + a^2*b*x^3)*log(b*x^3 + a) + 18*(b^3*x^9 +
 2*a*b^2*x^6 + a^2*b*x^3)*log(x))/(a^4*b^2*x^9 + 2*a^5*b*x^6 + a^6*x^3)

Sympy [F]

\[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=\int \frac {1}{x^{4} \left (\left (a + b x^{3}\right )^{2}\right )^{\frac {3}{2}}}\, dx \]

[In]

integrate(1/x**4/(b**2*x**6+2*a*b*x**3+a**2)**(3/2),x)

[Out]

Integral(1/(x**4*((a + b*x**3)**2)**(3/2)), x)

Maxima [A] (verification not implemented)

none

Time = 0.22 (sec) , antiderivative size = 117, normalized size of antiderivative = 0.62 \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=\frac {\left (-1\right )^{2 \, a b x^{3} + 2 \, a^{2}} b \log \left (\frac {2 \, a b x}{{\left | x \right |}} + \frac {2 \, a^{2}}{x^{2} {\left | x \right |}}\right )}{a^{4}} - \frac {b}{\sqrt {b^{2} x^{6} + 2 \, a b x^{3} + a^{2}} a^{3}} - \frac {1}{6 \, {\left (x^{3} + \frac {a}{b}\right )}^{2} a^{2} b} - \frac {1}{3 \, \sqrt {b^{2} x^{6} + 2 \, a b x^{3} + a^{2}} a^{2} x^{3}} \]

[In]

integrate(1/x^4/(b^2*x^6+2*a*b*x^3+a^2)^(3/2),x, algorithm="maxima")

[Out]

(-1)^(2*a*b*x^3 + 2*a^2)*b*log(2*a*b*x/abs(x) + 2*a^2/(x^2*abs(x)))/a^4 - b/(sqrt(b^2*x^6 + 2*a*b*x^3 + a^2)*a
^3) - 1/6/((x^3 + a/b)^2*a^2*b) - 1/3/(sqrt(b^2*x^6 + 2*a*b*x^3 + a^2)*a^2*x^3)

Giac [A] (verification not implemented)

none

Time = 0.30 (sec) , antiderivative size = 120, normalized size of antiderivative = 0.64 \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=\frac {b \log \left ({\left | b x^{3} + a \right |}\right )}{a^{4} \mathrm {sgn}\left (b x^{3} + a\right )} - \frac {3 \, b \log \left ({\left | x \right |}\right )}{a^{4} \mathrm {sgn}\left (b x^{3} + a\right )} - \frac {9 \, b^{3} x^{6} + 22 \, a b^{2} x^{3} + 14 \, a^{2} b}{6 \, {\left (b x^{3} + a\right )}^{2} a^{4} \mathrm {sgn}\left (b x^{3} + a\right )} + \frac {3 \, b x^{3} - a}{3 \, a^{4} x^{3} \mathrm {sgn}\left (b x^{3} + a\right )} \]

[In]

integrate(1/x^4/(b^2*x^6+2*a*b*x^3+a^2)^(3/2),x, algorithm="giac")

[Out]

b*log(abs(b*x^3 + a))/(a^4*sgn(b*x^3 + a)) - 3*b*log(abs(x))/(a^4*sgn(b*x^3 + a)) - 1/6*(9*b^3*x^6 + 22*a*b^2*
x^3 + 14*a^2*b)/((b*x^3 + a)^2*a^4*sgn(b*x^3 + a)) + 1/3*(3*b*x^3 - a)/(a^4*x^3*sgn(b*x^3 + a))

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{x^4 \left (a^2+2 a b x^3+b^2 x^6\right )^{3/2}} \, dx=\int \frac {1}{x^4\,{\left (a^2+2\,a\,b\,x^3+b^2\,x^6\right )}^{3/2}} \,d x \]

[In]

int(1/(x^4*(a^2 + b^2*x^6 + 2*a*b*x^3)^(3/2)),x)

[Out]

int(1/(x^4*(a^2 + b^2*x^6 + 2*a*b*x^3)^(3/2)), x)